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ABSTRACT 

EO/IR/Laser  detection  of  a  target  amidst  clutter/background  is  a  difficult  problem  often  treated  with  simplistic 
models.  Unlike  noise,  clutter  is  more  complex,  neither  spectrally  white  nor  statistically  Gaussian.  Therefore,  it  is 
insufficient  to  lump  clutter  with  noise  and  use  standard  detection  curves.  Using  current  target  detection  models,  it  is 
extremely  difficult  to  perform  effectiveness  assessments  of  signature  management  technologies  for  survivability  of  military 
ground  vehicles.  Current  models  do  not  consider  the  vehicle  on  a  component-level  basis  and  do  not  account  for  artifacts 
introduced  into  images  from  aliasing  and  varying  amounts  of  clutter.  Algorithms  must  be  developed  that  quantify  the 
effects  of  random  backgrounds  on  the  imaging  capability  of  electro-optical  systems  to  improve  false  alarm  rates.  Current 
trends  dictate  that  EO/IR/Laser  imaging  systems  must  consider  developments  in  signature  management  technologies  and 
countermeasures  that  are  driving  clutter  magnitudes  higher  than  target  signature  magnitudes.  These  trends  make  the 
problem  of  target  detection  in  clutter  especially  critical.  Battelle  has  produced  image  randomization  software  called 
BATRAN  (Background  and  Target  Randomization)  which  computes  various  types  of  statistical  distributions  to  randomize 
background  and  target  pixels  separately.  The  types  of  statistics  implemented  include  exponential,  Gaussian,  log-normal, 
and  Rice  distributions  for  both  the  background  and  target.  To  generate  synthetic  images  to  assess  the  detection 
performance  of  thermal  imaging  systems  and  countermeasured  platform  signatures,  a  method  to  characterize  the 
background  and  target  is  required  so  that  their  signatures  can  be  statistically  matched.  Current  methods  use  an  area- 
weighted  average  temperature  difference  (AWAA7),  which  is  regarded  as  inadequate  in  representing  observer’s  sensitivity 
to  the  inherent  detection  cues  of  the  target/background/clutter  signatures.  In  an  effort  to  identify  a  more  robust  and 
accurate  A T  metric  definition  for  background  and  target  matching,  Battelle  also  developed  a  new  A T  metric  definition  and 
its  equation  using  RMS  pixel-based  higher  order  statistics  for  the  background  and  target  signature  pixel  data  in  a  scene 
image.  This  new  AT  metric  provides  a  better  estimate  of  true  signature  difference  between  the  background/clutter  and 
target,  enabling  more  accurate  matching  of  the  background/clutter  and  target  for  use  in  sensor  detection  performance 
assessment. 


I.  INTRODUCTION 

BATRAN,  a  software  package  produced  by  Battelle,  randomizes  each  background  and  target  pixel  in  an  input  scene 
image.  The  randomizing  distribution  of  the  background  and  ta.get  pixels  can  be  different.  The  properties  of  randomizing 
statistical  distributions  are  provided  as  a  priori  distributions.  A  method  to  randomize  pixels  in  an  image  is  also  provided. 
Output  images  from  BATRAN  are  presented  along  with  inherent  input  images.  The  output  images  are  randomized  with 
various  pairs  of  randomizing  distributions.  These  output  images  are  then  used  for  AT  analysis. 

Currently  used  AT  metric  definitions  and  their  equations1 ,2,3  are  not  proper  representations  for  describing 
temperature  varying  features  of  a  target  and/or  its  background.  Some  AT  metrics  oversimplify  the  effects  of  these  features 
and  misrepresent  the  apparent  temperature  difference  which  provides  an  inherent  detection  cue  of 
target/background/clutter.  Other  AT  metrics  treat  the  pixels  from  a  target  and  its  background  in  a  sensor  image  as  if  they 
come  from  two  separate  images,  which  ignores  a  relative  temperature  apparentness  of  the  target  and  its  background. 
These  facts  inadequately  represent  a  human  observer’s  sensitivity  to  the  detection  of  target  from  background  and  clutter. 
Therefore,  these  problems  associated  with  the  currently  used  AT  metrics  necessitate  the  development  of  a  new  AT  metric 
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definition  and  its  equation.  In  this  paper,  a  new  A T  metric  is  defined.  An  equation  representing  this  new  definition  is  then 
derived.  This  new  AT  equation  is  implemented  in  BATRAN  along  with  other  AT  equations  with  which  the  new  AT  metric 
is  compared. 

The  remainder  of  this  paper  is  organized  as  follows:  Section  2  discusses  the  theory  of  the  background  and  target 
randomization  algorithm  implemented  in  BATRAN  including  the  a  priori  distributions  and  the  pixel  randomization  method, 
and  presents  the  randomized  output  images;  Section  3  discusses  the  new  AT  definition  and  its  equation  derivation,  and 
compares  the  new  AT  metric  with  other  AT  metrics;  Summary  and  conclusions  are  given  in  Section  4;  Section  5  contains 
acknowledgements;  and,  Section  6  lists  the  references  cited  during  this  research  performance  period. 

2.  BACKGROUND/CLUTTER  AND  TARGET  STATISTICS  MODELING  FOR  RANDOMIZATION 

This  section  discusses  the  theory  behind  the  background  and  target  pixel  randomization  algorithm  and  its 
development  and  modeling  in  BATRAN.  The  various  statistical  distributions  available  for  the  background  and  target 
randomization  are  provided  by  their  probability  density  function  (PDF)  and  cumulative  distribution  function  (CDF).  The 
random  number  generation  from  a  desired  statistical  distribution  using  the  inverse  transform  method  is  presented.  To 
assure  that  the  generated  random  numbers  represent  the  desired  statistics,  BATRAN-simulated  PDF  and  CDF  are 
compared  to  the  calculated  (or  exact)  PDF  and  CDF.  Finally,  the  BATRAN -randomized  output  images  are  displayed  to 
show  effects  of  the  randomization  by  different  statistics  on  a  couple  of  input  inherent  images. 

2.1  Properties  of  a  priori  distributions 

The  statistical  distributions  available  for  image  randomization  in  BATRAN  are  exponential,  Gaussian,  log-normal, 
and  Rice.  Because  these  distributions  are  not  derived  from  a  physical  description  of  clutter,  they  are  referred  to  as  a  priori 
distributions.  These  are  expressed  by  a  PDF  and  a  CDF,  denoted  by  p(x)  and  P(x),  respectively.  For  a  priori  distribution 
equations  presented,  p  and  o  denote  the  mean  and  standard  deviation  of  the  distribution,  respectively.  Another  parameter 
provided  for  the  log-normal  and  Rice  distributions  is  the  standard  deviation-to-mean  ratio,  denoted  by  R. 

2.1.1  Exponential  distribution 

The  PDF  and  CDF  of  the  exponential  distribution  are  given  in  (l.a)  and  (l.b),  respectively.  P(x)  is  the  definite 

X 

integral  of  p(x),  P(x)  =  jp{x)dx. 

p(x)  =  —  expf-— 1  (a);  P(x)  =  1  -  expf-— )  (h)  (1) 

P  V  Pj  V  pj 


Applying  (3)  to  (2),  (4)  gives  the  Gaussian  CDF  in  terms  of  the  error  function’s  CDF,  Pz(y)- 
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P(x)  = 


The  Gaussian  CDF  can  be  then  represented  as  in  (5),  which  is  the  error  function’s  CDF  approximation.  This  equation  is  a 
polynomial  and  rational  approximation4. 

P{x)  =  Pz(y,t)  =  1  -  Z(v)(hj/  +  b2t2  +  b3t3  +  b4t4  +  bsts )  +  e(v)  (5) 


where  v  -  ■■■—  *  and  t  -  -j— ^ — .  The  error  term  is  bounded  such  that  |  e(i>)  |  <7.5  X 1CT8.  The  coefficient  values  in  (5)  are 
bj=  0.319381530,  b2  =  -0.356563782,  b3  =  1.781477937,  b4  =  -1.821255978,  b5  =  1.330274429,  and  q  =  0.2316419. 

2.1.3  Log-normal  distribution5,4,7 


The  PDF  of  the  log-normal  distribution  is  given  in  (6). 


P(x)  = 


v/2^ouj 


exp 


flog.W-^taf 

—  - 


2  at 


(6) 


Since  the  log-normal  distribution  is  a  Gaussian  distribution  with  the  natural  logarithm  of  random  variables,  (5)  is 
used  with  v  =  '"l*  to  obtain  the  log-normal  CDF.  This  is  because  the  log-normal  CDF  can  be  written  as  (7). 


kv-”N 


\j2n 


/  exp 
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it 

>4 

°1- 

(7) 


Terms  in  (7)  are  <rln  =  logoff)  =  standard  deviation  of  log^j— j ,  and  =  \oge(m),  where  m  is  the  median  of  x.  The 
mean-to-median  ratio  (M  =  n /m)  of  x  is  related  to  the  standard  deviation  as  in  (8). 


ain  =  [2xlogf(Af)]J  =  2xlog,^J  1  =  [2x{logr(n)  -  m^}]1 


(8) 


The  median  can  then  be  found  by  solving  (8)  when  the  mean  and  the  standard  deviation  are  given  for  the  log-normal 
distribution  as  shown  in  (9). 


2 

™Tn  =  l°ge(m)  =  log»  -  ~ 


(9) 


The  standard  deviation-to-mean  ratio  (R)  is  expressed  in  (10. a),  which  can  be  solved  for  the  standard  deviation 
when  R  is  provided  as  indicated  in  (lO.b). 


CTia  =  /logj/?2  +  1)  (b) 


(10) 


2.1.4  Rice  distribution5,6 

The  PDF  of  the  Rice  distribution  is  expressed  in  (11),  where  I^a)  is  the  zero-order  modified  Bessel  function  of 
first  kind.  The  standard  deviation-to-mean  ratio  ( R )  is  expressed  in  (12. a).  When  R  is  given,  the  dominant-to-diffuse  ratio 
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(r)  can  be  expressed  as  in  (12. b).  The  dominant-to-diffuse  ratio  is  the  ratio  of  the  single  dominant  scatterer  signature  to 
the  total  signature  of  the  small  scatterers. 


P(x) 


1+r 

P 


exp 


- 1 

hUr)i}'\ 

r(l  +r)  — 

N  p; 

(ID 


1  +  r 


(a); 


r  =  ~g  +  Jgxig- 1)  ( b ); 


g  =  1  -R-1  (0  (12) 


Numerical  integration  of  the  Rice  PDF  is  necessary  to  obtain  its  CDF  because  there  is  neither  an  analytic  closed-form 
expression  nor  a  polynomial  approximation  of  the  Rice  CDF. 

2.1.5  Distribution  example 


Figure  1  shows  an  example  of  the  four  distributions  defined  above.  Shown  is  each  distribution’s  complement 

X 

cumulative  probability  (CCP),  1  -  ^p(x)dx.  Each  distribution  has  a  mean  of  1CT2.  The  second  parameter  of  each 
distribution  is  shown  in  Table  1. 


Figure  1  Example  of  cumulative  distributions 


Table  1.  The  value  of  the  mean  and  the  second 
parameter. 


Distribution 

type 

Mean 

Second  parameter 

Input  to 
BATRAN 

Exponential 

ft  =  IO'2 

N/A 

II 

© 

M 

Gaussian 

ft  =  io-2 

a  =  IO’5 

a  =  IO’3 

Log-normal 

ft  =  IO'2 
(m  =  1.96XKT1 
or  =  -6.2343) 

a  =  6.08 
(or  aln  =  1.805) 

R  =  5.0 

Rice 

ft  =  IO'2 

r  =  10.9782 

R  =  0.4 

2.2  Inverse  transform  method*  (ITM)  to  randomize  an  input  scene  image 

The  inverse  transform  method  is  applied  to  generate  random  deviates  (or  random  numbers)  to  randomize  the 
background  and  target  pixels  using  a  desired  distribution.  This  process  is  modeled  in  BATRAN.  Figure  2  is  used  to 
illustrate  the  ITM  implementation  in  BATRAN  with  the  following  steps. 

Step  1:  Generate  a  desired  CDF  either  analytically  or  numerically,  which  is  denoted  by  P(x)  and  represented  by  a 

solid  line  in  Figure  2. 

Step  2:  Generate  uniform  random  deviates  between  0  and  1,  which  is  denoted  by  y. 

Step  3:  With  P(x)  and  y,  one-to-one  matching  of  y  and  x  can  be  performed  using  P(x)  to  obtain  x,  which  is  a 

random  deviate  from  the  desired  distribution. 

Applying  ITM  to  an  input  (inherent)  scene  image,  pixels  are  randomized  with  a  desired  statistical  distribution.  It  is  known 
that  a  pixel  belongs  either  to  the  background  or  to  the  target  prior  to  the  randomization.  Therefore,  the  background  pixels 
can  be  randomized  by  one  statistical  distribution  and  the  target  pixels  can  be  randomized  by  another  distribution.  The 
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random  numbers,  which  may  represent 
fluctuations  in  temperature,  intensity,  or 
radiance,  are  generated  and  added  to  the 
inherent  background  or  target  pixel 
values. 


2.4  Comparison  of  exact  distribution 

and  the  BATRAN  generated 

distribution 

In  this  section,  random  variates 
generated  (or  simulated)  by  BATRAN  are 
examined  to  determine  whether  they 
closely  follow  the  selected  distribution  by 
comparing  them  to  the  exact  (or 
calculated)  distribution.  Figures  3,  4,  5, 
and  6  show  the  PDF  for  exponential, 

Gaussian,  log-normal,  and  Rice 
distributions,  respectively.  Figures  7,  8, 

9,  and  10  show  the  CDF  for  exponential, 

Gaussian,  log-normal,  and  Rice 
distributions,  respectively.  In  each  figure, 
the  solid  and  dotted  lines  represent  the 
calculated  and  simulated  distributions, 
respectively.  The  simulated  distribution  is 
generated  by  performing  105  trials.  From 
Figures  3,  4,  5,  and  6,  it  is  apparent  that 
the  number  of  trials  to  generate  the  simulated  distribution  is  106  since  there  are  no  samples  taken  beyond  the  PDF  of  ICT6. 
It  is  also  noted  that  as  the  simulated  PDF  approaches  the  tail  region  of  the  distribution,  it  deviates  from  the  calculated 
PDF.  This  is  because  not  enough  samples  are  taken  at  each  end  of  the  tail  regions  to  replicate  the  calculated  PDF. 
However,  the  simulated  CDF,  as  shown  in  Figures  7,  8,  9,  and  10,  agrees  with  the  calculated  CDF  up  to  the  complement 
cumulative  probability  (CCP)  of  ICT3.  Between  the  CCP  of  10"5  and  10-6,  the  simulated  CDF  deviates  from  the  calculated 
CDF.  Beyond  the  CCP  of  Kf6,  no  data  is  available.  Figures  3  and  7  show  the  exponential  PDF  and  CDF,  respectively. 
The  mean  (p)  of  this  distribution  is  10.  Figures  4  and  8  show  the  Gaussian  PDF  and  CDF,  respectively.  The  mean  (p)  is 
0  and  the  standard  deviation  (a)  is  20.  Figures  5  and  9  show  the  log-normal  PDF  and  CDF,  respectively.  The  mean  (p ) 
is  5  and  the  standard  deviation-to-mean  ratio  (R)  is  1.  From  the  mean  and  the  standard  deviation-to-mean  ratio,  the 
median  ( m )  and  the  standard  deviation  (<r)  are  3.5355  and  2.2992,  respectively,  using  (9)  and  (10. b).  In  natural  logarithm, 
the  median  and  the  standard  deviation  (o^)  are  1.2629  and  0.8326,  respectively.  Figures  6  and  10  show  the  Rice 
PDF  and  CDF,  respectively.  The  mean  ip)  of  the  distribution  is  35  and  the  standard  deviaiion-to-mean  ratio  (R)  is  0.5. 
From  the  standard  deviation-to-mean  ratio,  the  dominant-to-diffuse  ratio  (r)  is  6.4641  using  (12. b)  and  (12. c).  These 
figures  show  that  the  random  numbers  generated  in  BATRAN  closely  follow  the  exact  PDF  and  CDF.  This  not  only 
verifies  that  the  random  numbers  are  generated  with  the  input  statistics  parameters,  but  also  validates  that  they  are 
generated  by  the  desired  distribution. 

2.5  Output  images  from  BATRAN  after  background  and  target  pixel  randomization 

Example  output  images  from  BATRAN  after  randomization  were  generated  to  demonstrate  the  effects  of  the 
randomization  with  various  distributions.  A  simulated  and  a  measured  input  inherent  images  were  used  for  the 
randomization.  The  simulated  inherent  image  is  shown  in  Figure  11.  The  measured  inherent  image  is  shown  in  Figure 
12.  The  simulated  image  was  generated  by  the  Wright  Laboratories  Tactical  Decision  Aid  (TDA)  model  at  Battelle.  The 
measured  image  was  provided  by  TARDEC.  The  size  of  the  simulated  and  measured  images  were  378x377  and 
512x512  pixels,  respectively.  Background  pixel  values  of  the  simulated  inherent  image  were  set  to  zero;  therefore,  no 
difficulties  arose  in  separating  the  background  pixels  from  the  target  pixels.  However,  the  measured  image  had  to  be 


Figure  2  Illustration  of  inverse  transform  method  (ITM) 
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Probability  density  function  (PDF):  Calculated  vs.  Simulated 
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processed  to  separate  the  background  pixels  from  the  target  pixels.  The  Geographic  Resources  Analysis  Support  System 
(GRASS)  software  package  was  used  to  cut  out  the  target  pixels  from  the  background  pixels.  The  output  image  file  from 
GRASS  also  had  the  background  pixels  assigned  to  zero. 


Figure  11  Simulated  input  inherent  image.  Figure  12  Measured  input  inherent  image. 


2.5. 1  Output  image  examples  with  the  simulated/measured  image 

For  output  images  shown  in  Figures  13-16,  Table  2  summarizes  the  types  of  statistical  distributions  and  associated 
parameters  used  to  randomize  the  simulated  input  inherent  image.  For  output  images  shown  in  Figures  17-20,  Table  3 
summarizes  the  distributions  and  parameters  used  to  randomize  the  measured  input  inherent  image.  The  randomizing 
statistics  for  the  target  and  background  pixels  are  different  for  each  randomized  output  image.  These  figures  are  discussed 
in  detail  using  various  AT  metric  definitions  in  Section  3.3 

Table  2  Randomization  parameters  of  BATRAN  output  Table  3  Randomization  parameters  of  BATRAN  output 
image  using  the  simulated  input  inherent  image.  image  using  the  measured  input  inherent  image. 


Simulated 

Image 

Backgn 

Randomiz 

Dund  Pixel 
ing  Statistics 

Targ 

Randomiz 

et  Pixel 
ing  Statistics 

Figure  # 

Distribution 

Type 

Distribution 

Parameters 

Distribution 

Type 

Distribution 

Parameters 

13 

Exponential 

It  =  10 

Gaussian 

#t=0;  a=  10 

14 

Gaussian 

#i=24.45;cr  =  5 

Exponential 

It  =  5 

15 

Log-normal 

#t  =  6;  R  =  1 
mb  =  1.4452 
<jta  =  0.8326 

Rice 

#i=l;  R  =  0.5 
r  =  6.4641 

16 

Rice 

#t=15;  R=  1 

r  =  0 

Log-normal 

#1=12^=0.95 
mla  =  2. 1633 
<^  =  0.8020 

Measured 

Image 

Backgn 

Randomiz 

>und  Pixel 
ing  Statistics 

Targ 

Randomiz 

et  Pixel 
ing  Statistics 

Figure  # 

Distribution 

Type 

Distribution 

Parameter 

Distribution 

Type 

Distribution 

Parameter 

17 

Exponential 

##  =  25 

Log-normal 

It  =  7;  R  =  1 
=  1.5993 
=  0.8326 

18 

Gaussian 

#t  =  10 

o  =  10 

Rice 

##=  10;  /?=0.5 
r  =  6.4641 

19 

Log-normal 

#i=10;  /?= 0.9 
mta  =  2.0059 
0^  =  0.7703 

Gaussian 

H  =  7;  <j=10 

20 

Rice 

M=34J!=0.95 
r  =  0.45401 

Exponential 

#t  =  30 
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Figure  13  Exponential/Gaussian  randomization  for 
background/target  using  simulated  image. 


Figure  14  Gaussian/exponential  randomization  for 
background/target  using  simulated  image. 


Figure  15  Log-normal/Rice  randomization  for  Figure  16  Rice/log-normal  randomization  for 

background/target  using  simulated  image.  background/target  using  simulated  image. 
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Figure  17  Exponential/log-normal  randomization  for 
background/target  using  measured  image. 


Figure  18  Gaussian/Rice  randomization  for 

background/target  using  measured  image. 


Figure  19  Log-normal/Gaussian  randomization  for  Figure  20  Rice/exponential  randomization  for 

background/target  using  measured  image.  background/ target  using  measured  image. 
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3.  A  NEW  AT  DEFINITION 


In  this  section,  a  new  A T  metric  definition  is  stated,  and  the  equation  based  on  this  definition  is  derived.  This 
definition  is  compared  with  other  AT  metric  definitions  such  as  the  conventional  AT  (  =  AWA AT  ),  the  Night  Vision 
Laboratory  (NVL)  RMS  AT  in  [2],  and  the  AT  found  in  {3].  The  AT  results  from  these  various  definitions  are  also 
compared  by  using  Figures  13  -  20.  The  conventional  AT,  NVL  AT,  and  the  AT  from  [3]  are  denoted  by  ATconv,  ATiVVL, 
and  AT!3j,  respectively.  The  equations  for  A Tconv,  A T^,  and  AT{3J  are  given  in  (13. a),  (13. b),  and  (13. c),  respectively. 


~  Fr  Pa  ^ NVL 


jf  £  « 

iVT  i-l 


aT  (fc);  AT[3]  =  /(pr  -  pa)2  +  of  (c)  (13) 


where,  nT  :  The  mean  temperature  of  the  target  pixels 
Nt  :  The  number  of  pixels  on  target 
aT  :  The  standard  deviation  of  the  target  pixel  temperature 


\  Each  target  pixel’s  temperature 
HB  :  The  mean  temperature  of  the  background  pixels 


The  A Tconv,  as  shown  in  (13. a),  does  not  consider  temperature  varying  features  on  both  target  pixels  and 
background  pixels.  This  A  Tconv  definition  is  an  oversimplification  of  the  background  and  target  temperature  characteristics. 
The  A TjVVL  is  simply  a  standard  deviation  of  the  target  pixel  temperatures.  No  background  terms  are  involved  in  (13. b) 
indicating  the  background  temperature  is  ignored.  The  A T^j  considers  the  target  temperature  varying  feature  by  having 
Op  however,  the  background  temperature  varying  feature  is  not  considered.  It  is  evident  in  (13. b)  and  (13. c)  that  ATiWL 
and  ATpj  cannot  yield  a  negative  AT,  a  case  where  the  target  is  cooler  than  the  background.  These  deficiencies  indicate 
the  need  to  define  a  new  AT  metric  that  includes  higher  order  statistics  of  both  background  and  target  pixels’  temperatures. 


A  human  operator  does  not  see  a  separate  target  or  background,  but  rather  sees  both  at  the  same  time.  Therefore, 
an  overall  scene,  which  includes  the  target  and  background,  should  be  involved  in  defining  a  new  AT  metric.  This  is 
because  when  the  overall  scene  is  located  far  enough  away  such  that  the  entire  scene  falls  onto  a  pixel,  only  the  mean 
temperature  of  the  scene  is  realized.  However,  this  mean  temperature  is  produced  not  only  by  the  background  temperature 
but  also  by  the  target  temperature.  Therefore,  it  is  adequate  to  imitate  the  human  perception  of  the  scene  by  evaluating  the 
standard  deviation  of  the  target  pixels’  temperatures  conditioned  on  the  mean  temperature  of  the  scene  (=  0^5),  and  by 
evaluating  the  standard  deviation  of  the  background  pixels’  temperatures  conditioned  on  the  mean  temperature  of  the  scene 
(=  crB| 5).  Then,  a  new  AT  metric  is  defined  as  the  difference  between  the  standard  deviation  of  target  pixel 
temperatures  conditioned  on  the  mean  temperature  of  scene  and  the  standard  deviation  of  background  pixel 
temperatures  conditioned  on  the  mean  temperature  of  scene.  This  is  shown  in  (14). 


^  m.  ar|s  “  as]s 


(14) 


3.1  Formulation  of  the  new  AT  equation 

Prior  to  formulating  the  new  AT  metric  equation,  the  pixel  temperature  statistics  of  target/background  and  scene  are 
expressed  below.  (15. a)  and  (15. b)  show  the  mean  and  variance  of  the  target  pixel  temperatures,  respectively.  The  mean 
and  variance  of  pixel  temperatures  on  the  background  are  expressed  in  (16. a)  and  (16. b),  respectively. 


1  N* 

iyT  i-i 

(a); 

4-ij-E(p.r-^r 

N T  1 

-  i  bn1 

rVj. 

2 

-  Pr 

(b) 

(15) 

" B  i-l 

(a)-. 

i-i 

i-l 

!  2 

-  Pj 

(b) 

(16) 
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where,  PT  :  Pixel  temperature  on  target  P8  :  Pixel  temperatures  on  background 

o-f  :  The  variance  of  pixel  temperatures  on  target  aB  :  The  variance  of  pixel  temperatures  on  background 

HT  :  The  mean  of  pixel  temperatures  on  target  fiB  :  The  mean  of  pixel  temperatures  on  background 

Nt  :  Number  of  pixels  on  target  NB  :  Number  of  pixels  on  background 


The  mean  and  variance  of  the  overall  scene  are  given  in  (17. a)  and  (17. b),  respectively.  It  is  noted  that  Ns  =  NB  +  NT  . 


t= 1 


VtNt  +  »bNb 
Nt  +  Nb 


(a) 


a,  = 


f: 


z(r,s  -  4 .  ±  Erff  - 1  ^  - 14  -  ^  - 


(  nt  « , 

-  EpiB 


1. 1 


Nt  +  Nb 


(17) 


1-1 


i-1 


nt  +  nb 


(b) 


where,  Ps  :  Pixel  temperatures  on  scene 

Hs  :  The  mean  of  pixel  temperatures  on  scene 


a  2  :  The  variance  of  pixel  temperatures  on  scene 
Ns  :  Number  of  pixels  on  scene  (=  NB  +  NT) 


Then,  a ^s2  and  <rBiS2  can  be  expressed  in  terms  of  trp  N-p,  fiB,  aB,  NB,  and  Ns.  The  conditional  variances  are  shown 
in  (18. a)  and  (18. b)  for  the  target  and  background  pixel  temperatures,  respectively. 


Therefore, 
of  (18. b). 


4. 

nT  *  =  i 

4*  ■  i  S  ('* 


B  I»1 


the  new  AT  metric  equation  (  = 
This  is  written  in  (19). 


Ps  =  -\[HsX°2t  +  ~  Pa)2]  («) 

Ns  (18) 

Pr  =  -^*4  +  ^rx(Pa  "  Pr)2]  <*) 

Ns 

ATnrw)  is  the  difference  between  the  square  root  of  (18. a)  and  the  square  root 


<4  +  Pr  ~  2PrPj  + 

4  ♦  Pa  ~  2 + 


AT 


aT\s  ~  CTa|s 


Nr 


^N]o\  +  n\x[vt-  \i„f  -  pla\  *  Npx{\iB  -  nTf 


(19) 


In  the  above  equation,  it  is  apparent  that  the  temperature  varying  feature  on  the  target  and  background  pixels  is  included  in 
terms  of  their  variances.  Also,  the  number  of  pixels  on  target,  background,  and  scene  is  also  included,  which  shows  (19) 
is  a  pixel-based  temperature  difference  between  the  target  and  background.  Of  course,  the  mean  temperatures  of  the  target 
and  background  are  also  present  in  (19). 

3.2  Comparison  of  A Tnew  metric  to  other  AT  metric  equations 

The  other  AT  metric  equations  introduced  earlier  (A Tconv,  A Ty^,  and  A T^j)  will  now  be  compared  to  the  A T^ 
metric.  The  A Tconv  metric  is  evaluated  by  the  mean  temperatures  of  the  target  and  background  as  shown  in  (13. a).  It 
ignores  temperature  variations  of  the  target  and  background  pixels.  Thus,  the  variances  of  the  target  and  background  are 
zero,  G2  =  aB  =  0.  Under  these  conditions,  the  A Tnew  metric  reduces  to  in  the  A Tconv  metric  as  shown  in  (20). 

loj.,0,0,-0  “  TT  [77g x(Pr- Pa)  “  ^7- —  ~  TT  x(^J’  +  ^a) *(Pr“  P*)  =  Pr  ”  P«  ^ Tctmv  (20) 

"s  ‘'s 

The  A TlVVL  metric,  as  mentioned  before,  is  simply  a  standard  deviation  of  the  target  pixel  temperatures  as  shown  in  (13. b). 
It  does  not  consider  any  background  characteristics  nor  the  mean  target  temperature.  These  conditions  can  be  written  such 
as  ng  =  aB  =  Ng  =  nT  =  0.  Under  these  conditions,  the  A Tnew  metric  reduces  to  the  ATy^  metric  as  shown  in  (21). 


1  1 
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(21) 


AT 


yf-o,0j-° 

Ng*  0,Hy=0 


=  oT=  AT, 


m 


The  AT^y  metric,  from  the  equation  shown  in  (13.c),  considers  the  "temperature  variance  on  target  pixels  but  not  on 
background  pixels.  It  also  considers  the  mean  temperature  of  the  target  and  background  as  the  A Tconv  metric.  This 
condition  can  be  written  as  aB  =  0,  which  means  a  uniform  background.  For  this  condition,  the  AT^  metric  becomes 
(22). 

AT_|^0  =  -U/tfaVr  ~  ~  Ml  (22) 

"s 


It  seems  that  there  is  no  direct  relationship  between  (22)  and  the  A Tpj  metric  unless  some  special  conditions  are 
considered.  If  the  mean  of  the  target  and  background  pixel  temperatures  are  the  same  (p.T  =  nB),  (13. c)  and  (22)  become 
the  standard  deviation  of  the  pixel  temperatures  on  the  target  only,  which  is  the  ATyy^  metric  as  shown  in  (23). 


AT 


«,=0,t»T»|i,  "  A7[3lLr-|i,  "  A7WZ. 


=  0. 


(23) 


If  the  variance  of  the  target  pixel  temperatures  (ay2)  is  zero,  i.e.,  a  uniform  target  temperature,  (13. c)  and  (22)  become  the 
difference  between  the  target  mean  temperature  and  the  background  mean  temperature.  Thus,  in  this  case,  AT  is  the 
A  Tcanv  metric  as  shown  in  (24). 


AT. 


new  I  <3j  =  0,  oT=0 


=  AT, 


[31  1 


-  atcotv  -  pr  - 


(24) 


Therefore,  the  A  Tconv,  A T^,  and  A T[3]  metrics  describe  special  cases  of  the  ATn<w  metric.  The  A Tnew  metric  given  in 
(19)  includes  various  effects  due  to  the  target  and  background  characteristics  in  a  temperature  scene  image  by  evaluating  a 
number  of  pixels  on  the  target  and  the  background,  the  mean  of  the  target  and  background  temperatures,  and  the  variance 
of  the  target  and  background  temperatures.  The  variance  represents  the  temperature  varying  feature  of  the  target  and 
background.  The  ATnew  metric  is  derived  from  a  conditional  standard  deviation  of  the  target  and  background  pixel 
temperatures  on  an  overall  scene  as  shown  in  (18. a),  (18. b),  and  (19).  However,  the  other  AT  metrics  are  derived  by 
evaluating  the  target  and  the  background  separately  as  if  they  belong  to  separate  images.  This  is  not  a  correct 
interpretation  of  the  image  under  evaluation,  since  an  operator  views  the  target  pixels  and  background  pixels  in  the 
evaluated  image  simultaneously. 

3.3  A Tnew  metric  comparison  with  other  AT  metrics  using  Figures  13  through  20 

Figures  13  -  20  are  used  to  obtain  AT  values  using  the  A Tconv,  ATiVVL>  AT[3j,  and  A Tnew  metric  equations  given  by 
(13. a),  (13. b),  (13. c),  and  (19),  respectively.  Table  4  lists  the  target  and  background  statistics  and  ATs  of  Figures  13  - 
16.  Table  5  shows  the  same  for  Figures  17  -  20. 

Using  images  in  Figures  13  -  16,  the  various  AT  definitions  discussed  in  Sections  3.0,  3.1,  and  3.2  were  evaluated. 
All  of  the  target  and  the  background  pixels  contributed  to  the  AT  results.  Observing  Figure  13  and  its  AT  values,  it  is 
difficult  to  tell  which  AT  represents  Figure  13.  Nevertheless,  the  AT  above  10  degrees  Kelvin  seems  too  high  compared  to 
Figure  13.  The  A T^  seems  to  be  an  adequate  result  for  Figure  13.  In  Figure  14,  the  background  mean  temperature  and 
the  target  mean  temperature  were  intentionally  matched  using  BATRAN  to  evaluate  the  performance  of  each  AT  metric. 
As  seen  in  Figure  14,  much  of  the  target  features  are  buried  in  the  background  due  to  the  matching.  The  A Tconv  is  zero, 
which  is  not  a  true  statement  because  the  target  can  still  be  detected  and  recognized.  The  ATy^  and  A Tpj  may  be 
reasonable  only  if  hot  parts  of  the  target  are  considered.  Since  all  of  the  target  pixels  were  considered  in  evaluating  AT, 
the  results  from  these  definitions  may  not  properly  describe  the  image.  The  result  from  the  A T^  metric  provides  a  more 
correct  representation  for  the  image.  In  Figure  15,  not  only  were  the  background  mean  temperature  and  target  mean 
temperature  closely  matched,  but  also  the  background  variance  and  the  target  variance  were  closely  matched  using 
BATRAN.  As  seen  in  Figure  15,  the  target  image  is  deteriorating.  Except  for  a  few  hat  parts,  most  of  the  target  parts 
are  relatively  cooler  than  the  background.  The  results  from  the  ATiVVL  and  A Tpj  metrics  are  only  the  standard  deviation 
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(<r)  of  the  target  pixel  temperatures,  and  the  results  from  these  metrics  are  relatively  high  compared  to  the  image.  The 
A Tnm  and  ATconv  values  are  very  close  to  each  other.  This  is  because  the  variance  of  the  background  and  the  target  pixels 
were  matched.  These  results  may  not  reflect  target  recognition  ability.  In  Figure  15,  human  eyes  can  easily  detect  and 
recognize  the  target..  The  AT  results  using  the  image  depicted  in  Figure  16  show  that  the  A Tconv,  A and  A 7^  values 
are  too  high,  all  being  greater  than  12  degrees  Kelvin.  The  A value  shows  a  reasonable  result  for  the  image. 

Table  4  Simulated  image  output  pixel  statistics  and  Table  5  Measured  image  output  pixel  statistics  and 

various  ATs  comparison.  various  ATs  comparison. 
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Using  images  in  Figures  17  -  20,  the  various  AT  metrics  are  also  discussed.  The  background  mean  temperature 
and  the  target  mean  temperature  are  closely  matched  using  BATRAN  in  Figure  17.  This  results  in  A Tconv  close  to  zero. 
The  A TfjvL  and  AT^  values  are  the  standard  deviation  of  the  target  pixel  temperatures,  and  their  results  are  too  high 
compared  to  Figure  17.  The  overall  target  temperature  is  cooler  than  the  overall  background.  This  is  shown  in  ATnfw 
with  a  negative  value.  As  previously  noted  in  Section  3.0,  the  AT^  and  ATpj  metrics  cannot  discern  whether  the 
background  is  cooler  than  the  target,  or  vice  versa.  Figure  18  closely  represents  the  image  taken  during  30  mm/hour  rain 
fall.*  The  BATRAN  randomization  can  be  used  to  generate  images  which  may  represent  the  target  and  background  under 
various  weather  conditions.  The  AT  results  from  Figure  19  show  that  the  A 7]^  and  A T^j  values  are  too  high  compared 
to  the  image.  Either  A Tconv  or  A Tnew  value  represents  a  reasonable  temperature  difference  between  the  background  and  the 
target.  The  AT  results  calculated  from  Figure  20  clearly  tell  the  tolerance  of  each  AT  metric.  The  image  may  represent 
the  background  and  the  target  under  a  severe  weather  condition.  The  target  is  barely  distinguishable  in  the  image. 
However,  the  A Tconv,  A 7^,  and  AT[3j  metrics  yield  the  temperature  differences  of  13.6,  31.69,  34.48  degrees  Kelvin, 
respectively.  These  unreasonable  results  were  found  because  the  definition  of  these  metrics  considers  the  target  and  the 
background  separately  as  if  they  belong  to  separate  images.  Also,  these  definitions  do  not  consider  the  background 
temperature  varying  features  in  their  equations.  However,  the  ATnew  metric  interpreted  the  image  under  a  severe  weather 
condition,  as  shown  in  Figure  20,  quite  well. 

4.  SUMMARY  AND  CONCLUSIONS 

Background,  clutter,  and  target  field  measurements  from  EOAR,  ladar  (laser  radar),  visual,  and  other  sensor 
systems  show  their  statistics  can  be  represented  by  the  distributions  implemented  in  the  BATRAN  (Background  and  Target 
Randomization)  software,  namely  exponential,  Gaussian,  log-normal,  and  Rice  distributions.  The  background  and  target 
pixel  randomization  algorithm  using  these  distributions  was  first  developed  analytically  prior  to  implementation  in  the 
BATRAN  software,  which  performs  the  randomization  on  an  input  scene  image.  The  background  can  be  randomized  by  a 
statistical  distribution  which  differs  from  the  target  randomizing  statistical  distribution.  Also,  background  and  target 
randomizing  statistical  distributions  may  be  the  same.  The  BATRAN  software  allows  for  insertion  of  randomness,  which 
may  represent  artifacts  due  to  aliasing,  background  characteristics,  and  varying  amounts  of  clutter,  into  sensor  images. 
The  randomness  is  produced  according  to  the  given  statistical  distributions  so  that  the  randomized  image  can  be  used  to 


Figure  18  was  compared  to  a  real  image  taken  by  TARDEC  during  30  mm/hr  rain  fall. 
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analyze  various  detection  and  ATR  algorithms.  BATRAN  is  also  able  to  incorporate  weather  effects  on  sensor  image  if 
the  weather  statistical  characteristics  are  known.  With  these  capabilities,  BATRAN  is  a  tool  for  enhancing  the 
effectiveness  assessment  of  signature  management  technologies  for  military  ground  vehicles,  thus  increasing  their 
survivability. 

The  background  RMS  matching  with  a  target  was  performed  with  a  newly  developed  AT  metric  definition.  This 
new  AT  metric,  denoted  by  A T^  in  this  paper,  performed  better  than  the  other  definitions  reviewed  in  this  paper  when 
applied  to  various  scene  images,  both  simulated  and  measured.  Once  the  background/clutter  and  target  statistics  are 
known,  and  can  be  described  any  one  of  the  BATRAN  implemented  statistics,  the  BATRAN  software  can  match  the  target 
signature  to  the  background/clutter  well  by"  adjusting  the  given  target  and/or  background  statistical  parameters  to  provide  a 
desired  temperature  difference  (AT).  Various  combinations  of  background/clutter  and  target  statistics  generated  the  output 
images  (Figures  13  -  20)  from  the  input  scene  images  (Figures  11  and  12).  These  output  images  indicate  each 
combination’s  ability  to  match  a  target  to  the  background/clutter  to  a  human  observer.  Using  the  AT^  metric,  the  target 
and  background  can  be  matched  not  only  in  their  mean  temperatures  but  also  in  their  temperature  varying  features,  i.e., 
variance  or  standard  deviation.  This  corrects  the  deficiencies  of  the  existing  AT  metrics  mentioned  in  Section  3.0.  It  was 
also  concluded  that  the  A T^  metric  gave  the  best  apparent  match  of  the  suppressed  vehicle  with  the  background.9 
Therefore,  the  A T^  metric  allows  for  a  more  robust  and  accurate  assessment  of  the  detection  performance  of  thermal 
imaging  systems  than  the  other  AT  metrics. 
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